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Oh. Abstract 

C^ , Electromagnetic waves, in vacuum or dielectrics, can be confined in 

D ■ unbounded cylinders in such a way that they turn around the main axis. 

OlJ' For particular choices of the cylinder's section, interesting stationary con- 

CZ2 , figurations may be assumed. By refining some results obtained in previous 

O ■ papers, additional more complex situations are examined here. For such 

C/3 | peculiar guided waves an explicit expression is given in terms of Bessel's 

^~v functions. Possible applications are in the development of whispering 

i-Q ' gallery resonators. 
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00 ; 1 Introduction 

The classical equations of electromagnetism allow for solutions confined in ring- 
shaped domains. In vacuum, this is made possible by the orthogonality of 
the electric and magnetic fields (E and B) and by the enforcement of the two 
divergence-free conditions (divE = and divB = 0). In fact, if the lines of force 
are closed and orthogonal, a toroid is a natural environment to set up the initial 
conditions, that successively evolve according to: 



<9E , <9B 

H ■ -7T- = c 2 curlB — = - curlE 1) 

03 ■ dt at w 

where c denotes the speed of light. The above time-dependent equations are 
easily put in relation with the vector wave equation and its corresponding eigen- 
modes. 

The search for electromagnetic waves trapped in a toroid poses interesting 
mathematical questions. Numerical computations show a variety of solutions, 
whose dynamics depends on the section's shape. The behavior is strikingly sim- 
ilar to that of a non-viscous fluid confined in a vortex ring, but with additional 
intrinsic constraints. 



Explicit full solutions in terms of Bessel functions are available in the case of 
cylinders, where the magnetic field oscillates parallel to the axis and the electric 
field lays on the circular sections. The configuration recalls that of a train of 
solitons smoothly circulating inside a 2D rounded cavity. As documented in pQ , 
thin rings with large diameter and circular section can be well approximated by 
the above mentioned solutions. 

For more compact rings the use of numerical simulations is a necessity. By 
the way, not all the shapes are workable. Indeed, only a restricted range of 
sections are compatible with the electromagnetic constraints. Thus, the solution 
process must be implemented together with a sort of shape-detection algorithm 
(see[T]). 

Extensions of previous results (briefly recalled in the next section) are here 
obtained for electromagnetic waves trapped in unbounded cylinders where, the 
rotation around the axis is combined with a radial oscillation. Here we assume 
that the section is an annulus, so that the corresponding domain is a hollow 
cylinder. Having in mind the vector wave equation, the study is connected to 
the search of eigenfunctions of the Laplace operator in such a way that the 
dimension of the corresponding eigenspace is equal to four. As we will see, this 
analysis leads to the study of specific properties of Bessel's functions. 

2 Electromagnetic waves turning around an axis 

We assume that k > 1 is an integer number. We recall that the k-th Bessel's 
functions, of the first and the second kinds respectively, are independent solu- 
tions to the same eigenvalue problem (r > 0): 

J'^r) + ^ - ^M = _^ Jk{ur) (2) 

n"M + ^ - *^) = -, 2 nn ( 3) 

for a given parameter w. We also recall that Jk tends to zero for r — > 0, while 
Yfc is unbounded in the neighborhood of r = (see, e.g., [5]). 

Solutions of the entire set of Maxwell's equations are obtained in cylindrical 
coordinates (r, z,(f>) as follows (see [2], O): 

/ kf(ior) 
E = I cos(cu;£ — kip), 0, f'(ojr) sin(co;t — fc< 

B = -(o, f(ujr)cos(cut-k<P), o) (4) 

where / is a linear combination of Jk and Yfc. The magnetic field is parallel to 
the cylinder's axis and the electric field belongs to the plane (r, </>). Bounded 
solutions are possible in two cases: the domain is such that r > R m > for 
some given R m (see for instance the second picture in Fig. 1); the function / is 



just a multiple of Jk , thus r — may be included in the domain (see for instance 
the first picture in Fig. 1). 

The displacement of the electric field at a given time is shown in Fig. 1 for 
k = 1. During time evolution, the electromagnetic wave rigidly rotates around 
the central point. Two examples are taken into account. In the left picture, 
the domain is the circle < r < Rm, where Rm is the first zero of J\, which 
takes approximately the value 3.832. In this way, at the boundary, both the 
magnetic field and the radial component of the electric field are zero. In the 
right picture, we have the annular region R m = 1 < r < Rm- The outer 
diameter is now chosen in such a way that the electric field is orthogonal to the 
boundary (perfect conductivity). 
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Figure 1: Displacement of the electric field at fixed time obtained from the 
expression in |4]) subject to different kinds of boundary constraints. The sizes 
of the domains are not arbitrary, but strictly related to the frequency of rotation 
as well as boundary conditions. The magnetic held is orthogonal to the page. 



Once the frequency of rotation (depending on clj) is fixed and the type 
of boundary conditions decided, the size of the domain is automatically con- 
strained. This means that there are few domains allowing for the creation of 
these waves. If the proper size is not respected, the wave will not follow a sta- 
tionary motion, i.e., it does not complete the cycle preserving the phase. The 
conditions that permit such a guided behavior depend on the eigenvalues of the 
Laplace operator on the domain. If the shape is such that some eigenvalues 
have multiplicity higher than one, then the construction of such rotating soli- 
tary waves is possible (we show later the determination process). An analysis 
of spherical vortices has been provided in [I], exactly with the purpose of de- 
tecting domains where two independent eigenfunctions are related to the same 



eigenvalue. In this paper, we stay for simplicity in the case of the cylinder 
by analyzing the peculiar case where four independent eigenfunctions share the 
same eigenvalue. The evolution in a thin toroid can be approximated, with a 
rather good level of accuracy, starting from the cylinder's version; this extension 
however is not going to be studied here. 

3 A more involved evolution 

We would like to analyze the situation in which the multiplicity of one of the 
eigenvalues of the Laplace operator, on a certain annular domain, is equal to 
four. This will allow us to build more complex electromagnetic waves turning 
around an axis. We will need to play with both the functions Jk and Y^, 
therefore we have to stay away from the point r = (recall that Yjt is singular 
there). From now on we assume that R m = 1 < r < Rm, where the major 
diameter Rm has to be properly determined. Of course, if R m is modified the 
entire setting scales accordingly due to the linearity of Maxwell's equations. 

We impose Dirichlet boundary conditions on both boundaries (the inner and 
the outer circumferences), though other conditions may be considered. Taking 
Dirichlet conditions in the inner part involves working with the function: 

fn{r) = J„(Ar)Y„(A) - Y n {\r)J n {\) (5) 

where A > is a parameter. Of course, we have / n (l) = 0, Vn and VA. 
Now, we would like to have: 

jfn(R M ) =0 

< m * n (6) 

\fm(RM) = T W 

which is a nonlinear problem leading to the detection of Rm and A. This 
means that, once the diameter of the inner circumference has been fixed, the 
frequency of the evolving wave and the magnitude of the entire domain are going 
to be uniquely determined by (j6|) . Not necessarily such a problem has solution. 
However, with the help of numerical tests, we were able to establish some facts. 

Here below we report some of the conclusions of our analysis: 

m = 1, n = 2 : no values of Rm and A were found 

m = 1, n = 3 : no values of Rm and A were found 

m = 1, n = 4 : R M = 7.0927 for A = 1.0698 

m = 1, n = 5 : R M = 3.6761 for A = 2.3871 

to = 1, n = 6 : R M = 2.7603 for A = 3.6054 

m = 2, n = 3 : no values of Rm and A were found 

m = 2, n = 4 : no values of Rm and A were found 

to = 2, n = 5 : R M = 5.1500 for A = 1.7032 



Let us better examine a specific situation (m = 1, n = 4). According to Fig. 
2, both functions /i and fi are zero for r = R m = 1 and Rm ~ 7.0927. Thus, by 
imposing the same boundary constraints, we are able to find out two branches 
of vector solutions of the type given in @ (k = m and k = n respectively) 
related to the same frequency of rotation. These are generated by magnetic 
fields whose intensity is given by the level lines of Fig. 3. 




Figure 2: Plots of f\ and Ja for A 
vanish at r = 1 and r ss 7.0927. 



1.0699. In this case both the functions 



Successively, starting from A, we can build general solutions of the form: 

w\ sin(cvA t — <po) + W2 sin(cvA t) + w^ cos(cvA t — 0o) + 1/74 cos(cvA t) (7) 

The function Wj, j — 1, 2, 3, 4, are schematically reported in Fig. 4. In particu- 
lar, w\ and W3 are as shown in the first picture of Fig. 3, but with a difference 
of phase of 90 degrees. Combining w\ and W3 we can obtain the electromagnetic 
fields according to the expression (|4]) for k = 1. The other two functions, u>2 
and W4, are as shown in the second picture of Fig. 3 and differ for an angle of 
22.5 degrees. They also lead to ((4]) (this time with k = 4). The phase lag 4>o is 
arbitrarily given. 

It is just a direct computation verifying that the expression provided in ([7]) 
solves the wave equation in vacuum. In fact, the second derivative in time 
produces the multiplicative factor — c 2 A, while the application of the Laplace 
operator is equivalent to a multiplication by — A = — us 2 (see ([2]) and ([3])). 



In Fig. 5, one can see, at different time steps, the evolution of (|7|) for O = 0. 
Only half cycle is displayed. From the last picture the sequence restarts from 
the beginning with the color interchanged. An animation is available in [3] . The 
effect vaguely recalls the juggling of three clubs. 





Figure 3: Plots of two independent eigenfunctions of the Laplace operator with 
homogeneous Dirichlet boundary conditions on an annular domain. Despite the 
fact that they look totally different, the domain has been dimensioned in such 
a way the corresponding eigenvalues are the same. 
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Figure 4: Schematic representation of four independent eigenfunctions on the 
annular domain having the same eigenvalue. 
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Figure 5: Evolution at different times of a linear combination of eigenfunctions 
corresponding to the same eigenvalue. The two eigenfunctions of Fig. 3 merge 
and rotate with a periodic behavior. The frequency is uniquely established by 
the global size of the domain. 



4 Conclusions 

From the theoretical viewpoint, Maxwell's equations allow for special solutions 
trapped inside an infinite cylinder or, more or less equivalently, in a toroid. The 
study here developed in vacuum reveals original configurations and may help for 
instance in the analysis of confined plasma (see, e.g., [5]). Other applications 
may be found in the field of the so called whispering gallery resonators. Waves 
trapped in these cavities are smoothly guided to circulate around by continuous 
reflection returning at the origin with the initial phase. Spherical, cylindrical 
and ring-shaped whispering galleries are commonly produced for a broad range 
of industrial applications. Typical areas of interest are in fiber telecommunica- 
tions or biosensing. The literature is very rich. We just mention a couple of non 
extremely specialized publications: [6], [7]. 
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